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In this paper we show that dark matter in the form of dense matter/antimatter nuggets could 
provide a natural and unified explanation for several distinct bands of diffuse radiation from the core 
of the Galaxy spanning over 12 orders of magnitude in frequency. We fix all of the phenomenological 
properties of this model by matching to x-ray observations in the keV band, and then calculate the 
unambiguously predicted thermal emission in the microwave band, at frequencies smaller by 10 
orders of magnitude. Remarkably, the intensity and spectrum of the emitted thermal radiation 
are consistent with — and could entirely explain — the so-called "wmap haze": a diffuse microwave 
excess observed from the core of our Galaxy by the Wilkinson Microwave Anisotropy Probe (wmap). 
This provides another strong constraint of our proposal, and a remarkable nontrivial validation. If 
correct, our proposal identifies the nature of the dark matter, explains baryogenesis, and provides a 
means to directly probe the matter distribution in our Galaxy by analyzing several different types 
of diffuse emissions. 
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I. INTRODUCTION 



In this paper we discuss a testable and well-constrained 
model for dark matter [1, 2, 3, 4]. In particular, we 
explain how microwave emissions are an inevitable con- 
sequence of our proposal, and test the model against a 
recent analysis of the Wilkinson Microwave Anisotropy 
Probe (wmap) observations that suggest an anomalous 
emission from the core of our Galaxy (dubbed the "wmap 
haze") [5, 6, 7, 8, 9, 10]. Despite having no free 
parameters — the model is completely fixed by observa- 
tions at scales some 10 orders higher — our proposal is 
consistent with these observations, and could even ex- 
plain the anomaly if it survives further scrutiny. This 
provides a highly nontrivial test of our proposal, which 
remains consistent with all known constraints. 

To calculate the emissions, however, requires a care- 
ful analysis of several fields of physics, and it is easy to 
lose track of the overall structure of the calculation. We 
organize the paper as follows: In section II, we present 
a short review of the dark-matter proposal, emphasizing 
the assumptions that underlie the model, the observa- 
tional constraints that the model must satisfy, and outlin- 
ing the mechanism by which the dark matter will radiate, 
thus rendering it observable. In section III, we summa- 
rize the calculation of the thermal microwave emission 
from the dark matter, and show in section IV how this 
is consistent with the observations. Finally, in section V, 
we review all of the observational constraints of our pro- 
posal and discuss the testable predictions that it makes. 
To keep the logic clear, some technical details have been 
omitted from core of the paper. We include these in the 
appendix, completing our calculation, clcarpagc 
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II. PROPOSAL 
A. Dark Matter as Dense Quark Nuggets 

Two of the outstanding cosmological mysteries — the 
natures of dark matter and baryogenesis — might be ex- 
plained by the idea that dark matter consists of Com- 
pact Composite Objects (ccos) [1, 2, 3, 4] similar to 
Witten's strangelets [11]. The basic idea is that these 
CCOs — nuggets of dense matter and antimatter — form at 
the same QCD phase transition as conventional baryons 
(neutrons and protons), providing a natural explanation 
for the similar scales f2 DM as 50 B . Baryogenesis proceeds 
through a charge separation mechanism: both matter 
and antimatter nuggets form, but the natural CP viola- 
tion of the so-called 8 term in QCD 1 — which was of or- 
der unity 9 ~ 1 during the QCD phase transition — drives 
the formation of more antimatter nuggets than matter 
nuggets, resulting in the leftover baryonic matter that 
forms visible matter today (see [2] for details). Note, it 
is crucial for our mechanism that CP violation can drive 
charge separation. This idea may already have found ex- 
perimental support through the Relativistic Heavy Ion 
Collider (rhic) at Brookhaven [16]. 

This mechanism requires no fundamental baryon 
asymmetry to explain the observed matter/ antimatter 
asymmetry. From this, and the observed relation i7 DM w 
5fi B (see [17] for a review) we have 

^universe — -^nugget ~t~ -^visible -^antinugget (^"^) 
-^dark-matter — -^nugget ~t~ -^antinuggct ~ 5-Z? v isiblc (-^) 

where -B un ivcrse is the total number of baryons minus the 
number of antibaryons in the universe, -Bdark-matter is the 
total number of baryons plus the total number of an- 
tibaryons hidden in the nuggets and antinuggcts that 
make up the dark matter, S n uggot is the total number 
of baryons contained in all of the dark-matter nuggets, 
-Bantinugget is the total number of antibaryons contained 
in all of the dark-antimatter nuggets, and -B v isibie is the 
total number of residual "visible" baryons (regular mat- 
ter). Solving equations (1) gives the approximate ratios 

-Dantinugget--£?nugget*-£?visiblc — 3.2:1. 

Unlike conventional dark matter candidates, dark- 
matter/antimatter nuggets will be strongly interacting, 



1 If 8 is nonzero, one must confront the so-called strong CP prob- 
lem whereby some mechanism must be found to make the ef- 
fective 9 parameter extremely small today in accordance with 
measurements. This problem remains one of the most outstand- 
ing puzzles of the Standard Model, and one of the most natural 
resolutions is to introduce an axion field. (See the original pa- 
pers [12], [13], [14], and recent reviews [15].) Axion domain walls 
associated with this field (or ultimately, whatever mechanism re- 
solves the strong CP problem) play an important role in forming 
these nuggets, and may play in important role in their ultimate 
stability. See [1, 2] for details. 



but macroscopically large objects. They do not contra- 
dict any of the many known observational constraints on 
dark matter or antimatter [3, 18] for three reasons: 

(1) They carry a huge (anti)baryon charge \B\ as 10 20 
- 10 33 , so they have an extremely tiny number den- 
sity. This explains why they have not been di- 
rectly observed on earth. The local number den- 
sity of dark-matter particles with these masses is 
small enough that interactions with detectors are 
exceedingly rare and fall within all known detector 
and seismic constraints [3]. (See also [19, 20] and 
references therein. 2 ) 

(2) The nuggets have nuclear densities, so their in- 
teraction cross-section is small a/M as 10 -13 - 
10~ 9 cm 2 /g. This is well below typical astrophysi- 
cal and cosmological limits which are on the order 
of a/M < 1 cm 2 /g. Dark-matter-dark-matter in- 
teractions between these nuggets are thus negligi- 
ble. 

(3) They have a large binding energy such that baryons 
in the nuggets are not available to participate in 
big bang nucleosynthesis (bbn) at T as 1 MeV. In 
particular, we suspect that the baryons in these 
nuggets form a superfluid with a gap of the or- 
der A as 100 MeV, and critical temperature T c ~ 

as 60 MeV, as this scale provides a natural 
explanation for the observed photon to baryon ra- 
tio n B /n y ~ 10~ 10 [2], which requires a formation 
temperature of T form = 41 MeV [21]. 3 

Thus, on large scales, the nuggets are sufficiently dilute 
that they behave as standard collisionless cold dark mat- 
ter (ccdm). When the number densities of both dark 
and visible matter become sufficiently high, however, 
dark-antimatter-visible-matter collisions may release sig- 
nificant radiation and energy. In particular, antimat- 
ter nuggets provide a site at which interstellar baryonic 
matter — mostly protons and electrons — can annihilate, 
producing emissions with calculable spectra and energies 
that should be observable from the core of our Galaxy. 
These emissions are not only consistent with current ob- 
servations, but seem to naturally explain several mys- 
terious diffuse emissions observed from the core of our 
Galaxy, with frequencies ranging over some 12 orders of 
magnitude. 

Although somewhat unconventional, this idea natu- 



2 It is estimated in [19] that nuggets of mass from ~ 10 kg to 1 
ton (corresponding to B ~ 10 28 " 30 ) must account for less than an 
order of magnitude of the local dark matter. While our prefer- 
able range of B ~ in, 25-27 is somewhat smaller [18] and does not 
contradict [19], we still believe that B > 10 28 is not completely 
excluded by Apollo data, as the corresponding constraint is based 
on specific model dependent assumptions about the nugget mass- 
distribution [19], whereas nugget formation due to charge separa- 
tion as suggested in [2] may lead to a very different distribution. 

3 At temperatures below the gap, incident baryons with energies 
below the gap would Andreev reflect rather than become incor- 
porated into the nugget. 
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rally explains several coincidences, is consistent with 
all known cosmological constraints, and makes definite, 
testable predictions. Furthermore, this idea is almost en- 
tirely rooted in conventional and well-established physics. 
In particular, there are no "free parameters" that can 
be — or need to be — "tuned" to explain observations: In 
principle, everything is calculable from well-established 
properties of QCD and QED. In practice, fully calculat- 
ing the properties of these nuggets requires solving the 
fcrmion many-body problem at strong coupling, so we 
must resort to "fitting" a handful of phenomenological 
parameters from observations. 

Nevertheless, these unknown parameters may be de- 
termined to within an order of magnitude by observa- 
tions of processes at the keV scale (described below and 
in [4]). The model then makes unambiguous predic- 
tions about other processes ranging over more than 10 
orders of magnitude in scale. The point of this paper 
is to show that, remarkably, these unambiguous predic- 
tions are completely consistent with current observations, 
providing compelling evidence for our proposal, and ex- 
plaining another astrophysical puzzle: the origin of the 
so-called "wmap haze" [5, 6, 7, 8, 9, 10] . 4 We have con- 
sidered five independent observations of diffuse radiation 
from the core of our Galaxy: 

(1) Spi/integral observes 511 keV photons from 
positronium decay that is difficult to explain with 
conventional astrophysical positron sources [23, 24, 
25] . Dark-antimatter nuggets would provide an un- 
limited source of positrons as suggested in [26, 27]. 

(2) Comptel detects a puzzling excess of 1-20 McV 
7-ray radiation. We shall not discuss this here, but 
it has been shown in [28] that the direct c + c~ anni- 
hilation spectrum could nicely explain this deficit. 

(3) Chandra observes a diffuse keV x-ray emission 
that greatly exceeds the energy from identified 
sources [29]. Visible-matter/dark- antimatter anni- 
hilation would provide this energy. 

(4) Egret/crgo detects McV to GcV gamma-rays, 
constraining antimatter annihilation rates. We 
shall not discuss these constraints here, but it was 
shown in [4] that these constraints are consistent 
with the rates inferred from the other emissions. 

(5) Wmap has detected an excess of GHz microwave 
radiation — dubbed the "wmap haze" — from the in- 
ner 20° core of our Galaxy [5, 6, 7, 8, 9]. Annihila- 



4 We should remark here that our explanation of the WMAP haze 
with dark matter is not a new idea. It was suggested previ- 
ously that self-annihilating weakly interacting massive particles 
(WIMPs) might explain the WMAP haze [7, 8, 10]. These WIMPs 
must be very heavy, m ~ 100 GcV, and therefore, annihilation 
must produce significant amounts of high energy radiation [22] if 
it is to also explain microwave emissions with a typical frequency 
of u) ~ 10~ 4 . For example, if one takes central values for the 
WIMP parameters, then the microwave intensity from WIMP 
annihilation would be well below the observed intensity [22]. In 
any case, this proposal has very different predictions than ours. 



tion energy not immediately released by the above 
mechanisms will thermalize, and subsequently be 
released as thermal bremsstrahlung emission at the 
eV scale. Although the eV scale emission will be ob- 
scured by other astrophysical sources, the tail of the 
emission spectrum is very hard, and carries enough 
energy in the microwave to explain the wmap haze. 

To proceed, we start with the three basic postulates as- 
sumed in [4]: 

(A.l) The antimatter nuggets provide a virtually unlim- 
ited source of positrons (e + ) such that impinging 
electrons (e _ ) will readily annihilate at their sur- 
face through the formation of positronium [26, 27]. 

About a quarter of the positronium annihilations release 
back-to-back 511 keV photons. On average one of these 
photons will be absorbed by the nugget while the other 
will be released. 

(A. 2) The nuggets provide a significant source of an- 
tibaryonic matter such that impinging protons will 
annihilate. We assume that the proton annihilation 
rate is directly related to that of electrons through 
a suppression factor / < 1 as discussed in [4]. 

Proton annihilation events will release about 2 m p ss 
2 GcV of energy per event and will occur close to the 
surface of the nugget creating a hot spot that will ra- 
diate x-ray photons with keV energies containing some 
fraction g of the total annihilation energy. The remain- 
ing fraction 1 — g will be released at the eV scale, after 
the energy has thermalized within the nuggets. The tail 
of this thermal emission will be released in the microwave 
spectrum and may explain the observed WMAP haze. To 
further test this theory and connect all of these emissions, 
we make an additional assumption: 

(A. 3) We assume that the emitted 511 keV photons dom- 
inate the observed 511 keV flux, that the emit- 
ted keV x-rays dominate the observed diffuse x-ray 
flux, and that the thermally emitted microwaves 
dominate the observed WMAP haze. 

The basis for this assumption is that none of these fluxes 
has a convincing explanation. The nuggets may thus pro- 
vide the missing explanation in each case. This assump- 
tion allows us to use the observations at the keV scale — 
the 511 keV emission measured by integral, and the 
diffuse x-ray emission measured by Chandra — to fix all 
of the phenomenological parameters [4] . The model then 
makes unambiguous predictions about the properties of 
the WMAP haze, allowing it to be tested. 

As we shall see, the agreement is remarkable: even 
though our estimates are only accurate up to the order of 
magnitude, the picture that dark matter consists mostly 
of antimatter nuggets can completely explain all three of 
these puzzling emissions without contradicting a single 
observation. 
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B. Observable Dark Matter: Emissions 

As discussed in [4] , our proposal is that both electrons 
and protons annihilate on antimatter nuggets, releasing 
observable radiation from "hot spots" near the annihila- 
tion sites. These emissions are best though of as "jets" , 
and occur sufficiently rapidly that they are produced on 
a per-event basis, thus producing a spectrum that is in- 
dependent of the local environment [4] . This includes the 
511 keV spectrum from positronium annihilation [26, 27], 
the spectrum up to 20 MeV from direct c + c~ annihila- 
tion [28], the diffuse ~ 10 keV radiation from p + p~ an- 
nihilation [4], and the occasional GeV photon produced 
directly from proton annihilation [4]. 

The rates of annihilation and the energies released have 
been correlated [4] with the observed diffuse 511 keV 
positronium emission [23, 30, 31, 32, 33, 34] and dif- 
fuse ~ 10 keV emissions [29] observed from the core of 
the Galaxy, providing a test of the model. It was shown 
that both of these emissions could be nicely accounted 
for if the rate of x-ray energy released from p + p _ anni- 
hilation was related to that of c _ c + annihilation through 
a suppression factor / ■ j ~ 6 x 10~ 3 where the factor 
/ accounts for proton reflection from the sharp nuclear 
matter interface, and the factor ^ < g < 1 accounts 
for the fraction of the 2 GeV p + p~ annihilation energy 
released at the "hot-spots" . 

The topic of this paper is the remaining fraction 1 — g 
of the 2 GeV annihilation energy that will be transmit- 
ted deep withing the nuggets, ultimately being thermally 
radiated at a much lower energy scale. We shall show 
that most of this energy will be released at the eV scale, 
making it difficult to observe against the bright stellar 
background. The spectrum of this emission, however, 
will be shown to be extremely hard, resulting in a sig- 
nificant release of detectable microwave energy (~ 10~ 4 
eV). 

Our main point is that this microwave emission could 
fully account for the recently observed WMAP haze [5, 6, 
7, 35]: a puzzling diffuse emission from the core of our 
Galaxy. 

In Section III A we estimate the thermal emissivity of 
nuggets, and the spectrum of the emitted radiation. In 
Section III B we show how the nuggets reach thermody- 
namic equilibrium at an eV scale by balancing the anni- 
hilation rate with the emission. Armed with these esti- 
mates, in Section IV we compare the predictions of our 
model with the observations of the WMAP haze, using our 
previous results [4] to provide the normalizations, and ar- 
rive at the remarkable conclusion that our proposal natu- 
rally explains the energy budget IV A and spectrum IV B 
of the observations, even though the predictions are at an 
energy scale some 10 orders of magnitude smaller than 
the 511 keV scale at which the normalization was fixed! 
Finally, in Section VB we reiterate the testable predic- 
tions our model makes, thus providing a method with 
which to confirm or rule out the proposal over the next 
few years. 



III. THERMAL PROPERTIES OF THE 
NUGGETS 

A. Emissivity 

Here we discuss the properties of thermal emission 
from the "electrosphere" of the nuggets at low temper- 
atures T ~ eV. As we shall show in Section IIIB, this 
temperature can be established by comparing the rate 
of annihilation energy deposited in the nugget with the 
rate of emission. In what follows, we shall present a sim- 
ple estimate to capture the order of magnitude of the 
process. In principle, the exact numerical factors can be 
computed, but such a calculation is extremely tedious, 
and would be of no use since there are other uncertain- 
ties in this problem of a similar magnitude. 

The emissivity depends on the density n(z) of the 
positron cloud, which varies as a function of the dis- 
tance z from the quark matter core. At the eV scale 
temperatures, the most important region of emission will 
be the region of the electrosphere where the kinetic en- 
ergy of the particles p 2 /2m ~ T is on the same order 
as the temperature. Closer to the core of the nugget, a 
well-defined Fermi surface develops with pp T , and 
the low-energy excitations that can scatter and radiate 
arc confined to an effectively two-dimensional region of 
momentum space about the Fermi surface. As a result, 
there is a kinematic suppression of the emissivity from 
these regions and the emission will not change the order 
of magnitude estimate we present here. Sufficiently deep 
into the nugget, the plasma frequency will also be large 
enough that the emitted eV scale photons will be highly 
virtual and thus rapidly reabsorbed. A detailed discus- 
sion of the suppression of emission from the highly dense 
regions is presented in appendix 4. 

Here we shall estimate the emissivity of a Boltzmann 
gas of positrons. The Boltzmann approximation is valid 
where n <C p -3 ~ (mT) 3 / 2 and we can neglect both 
the fermion degeneracy that will suppress the emissivity 
closer to the core and many-body effects. Thus we start 
from the following expression 5 for the cross section for 
two positrons emitting a photon with to <C p 2 / (2m) [36], 

d^± a (^(l7+121n^)^. (2) 

15 \mJ \ mu> ) to 

The emissivity Q = dE/dt/dV — defined as the total en- 
ergy emitted per unit volume, per unit time — and the 



J Expression (2) should be contrasted with the well known dipole 
type of expression for different types of particles emitting soft 
photons, such as electrons and ions. With identical particles hav- 
ing the same charge to mass ratio e/m, the dipole contribution 
is zero, and the cross section is dominated by the quadrupolc 
interaction. This quadrupole character explains the appearance 
of the velocity (v) in numerator of (3) as opposed to the factor 
(v)~ 1 that enters the corresponding expression for electron-ion 
collisions. 
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spectral properties can be calculated from 

^-(u,z) = ni(z,T)n 2 {z,T)uj ( v 12 ^- 
auj \ auj 



4a / a \ 2 
15 



17+121n^i2. ) 

TOW / ' 



where n(z, T) is the local density at distance z from 
the nugget's surface, and v\2 = \v\ — V2\ is the relative 
velocity. The velocity and momentum p\2 need to be 
thermally averaged. To estimate this, we use the Boltz- 
mann ensemble at temperature T with a kinematic cutoff 
p 2 j (2m) > u>, as only particles with sufficient energy can 
emit photons with energy u>. In principle one can do bet- 
ter, but the current approach suffices to give the correct 
order of magnitude (sec Appendix 2 for details): 



Via 17 +12 In- 



J/T h(^) (4) 



where (this approximation is accurate to about 25%) 



h(x) 



17-121n(x/2) x < 1, 
17+121n(2) x>l. 



(5) 



To proceed with our estimates, we need the positron den- 
sity in the nugget's electrosphere at temperature T. As 
shown in Appendix 1, the corresponding expression in 
the nonrclativistic mean-field approximation is given by 



i(z,T) 



1 



T 



2ira (z + z) 2 ' 



(6) 



where z is a constant of integration to be determined by 
some appropriate boundary condition. It is known that 
the mean-field approximation is not valid for extremely 
large z, where exponential rather than power-law (6) de- 
cay is expected. We could accommodate the correspond- 
ing feature by introducing a cutoff at sufficiently large 
z = z ma x- The result, however, is not sensitive to this 
cutoff, so we shall simply take z max = oo below to obtain 
our order of magnitude estimate. 

Note: the electrosphere extends well beyond the core 
of the nugget. To see this, note that the Boltzmann 
regime (6) is based on an approximation that neglects 
the curvature of the nuggets surface (see appendix 1). 
This regime terminates only once the curvature becomes 
significant, i.e. once the electrosphere extends at least to 
the same order as the macroscopic radius of the nugget. 
Thus, the electrosphere occupies a significant fraction of 
the nugget's volume: it is not just a thin outer shell. Cor- 
rections from the finite size of the nugget are discussed 
in appendix 3 but do not affect the magnitude of our 
calculations. 

The parameter z is not a free parameter, but is fixed 
by matching the full density profile to the boundary of 



the nuclear matter core of the nugget, where the lepton 
chemical potential /i ~ 10 MeV is established by beta- 
equilibrium in the nuclear matter. A proper computation 
of z thus requires tracking the density through many or- 
ders of magnitude from the ultrarelativistic down to the 
nonrelativistic regime, which is beyond the scope of this 
work. The order of magnitude, however, is easily esti- 
mated by taking z = as the onset of the Boltzmann 
regime: 



T 



1 



n z=0 



27ra z 2 



(toT) 



3/2 



(7a) 
(7b) 



Numerically, z ~ 0.5 • 10 8 cm while the density n ~ 
0.3- 10 23 cur 3 for T ~ 1 eV. (See Appendix 1 for details 
of this calculation.) 

Our next task is to estimate the surface emissivity (ra- 
diant exitance) F = J dz Q(z) — defined as the energy 
E emitted per unit time dt, per unit area dA (flux) — 
from the nugget's surface by integrating the emissivity 
(3) over the Boltzmann regime z € [0, z max — * oo], and 
introducing an extra factor 1/2 to account for the fact 
that only the photons emitted away from the core can 
actually leave the system. 

Our final estimate for spectral surface emissivity can 
be expressed as follows: 



dF 



dE 



1 



^ (w) = dTdZd:"2X dz S (w ' z) 



4 T 3 a 5 / 2 4 /T 

45 7T 



1 



-u/T h ' 



T 



T 



(8) 



Integrating over u contributes a factor of T J dx (1 + 
x) exp(—x)h(x) ks 60 T, giving the total surface emissiv- 
ity: 



dE r dF 
/'tot = - . . = duj — — (a;) 



dt dA 



duj 



16 T 4 a 5/2 ,!~f~ 
— — y— (9) 

3 7T V TO 



Although a discussion of black-body radiation is inap- 
propriate for these nuggets (for one thing, they are too 
small to establish thermal equilibrium with low-energy 
photons), it is still instructive to compare the form of 
this surface emissivity with that of black-body radiation 
Fnn = aT 4 : 



F u 



320 



BB 



,5/2 



(10) 



At T = 1 cV, the emissivity F tot ~ 1Q~ 6 Fbb is much 
smaller than that for black-body radiation. The spectral 
properties of these two emissions are also very different 
at low frequencies u> < T as follows from (8). 

These two differences are essential to explain the WMAP 
haze. First, the suppressed total radiant exitance is re- 
quired to establish the eV temperature scale (this will be 
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discussed in Section IIIB). Second, the extremely long 
low-frequency tail due to the logarithmic dependence of 
h(x) is required to ensure that sufficient power is radi- 
ated in the microwave. Thus, as we shall show in IV, 
it is highly nontrivial that the scale of the emitted mi- 
crowave emission should be consistent with the observed 
WMAP haze emission: If the nuggets had been simple 
black-body emitters, the emission would be many orders 
of magnitude below the observed scale. 

Finally, we emphasize here that there are no free pa- 
rameters in this calculation: all of the scales are set 
by well-established nuclear and electromagnetic physics. 
The only unknown parameters that enter are the overall 
normalizations which can be fixed by considering the re- 
lated diffuse x-ray emission. This is the point of the next 
section. 



density in the bulge at a distance r ~ kpc from the core, 
where the observed WMAP haze originates: 



.5 kpc 



50. 



Combining these, we obtain 
10 9 GcV ( J_ 

W 



1 ) ' (l0" 3 c) 



( n \ 



V300/W 



(13) 



which must be compared with the total surface emissiv- 
ity (9) 



GcV / T 



F tot ~ 10 9 ^^- I — 1 



cm 2 • s V eV / 



\ 4+1/4 



B. Thermodynamic Equilibrium 

Armed with an estimate of the total emissivity (9), 
we may discuss the thermodynamic properties of the 
nuggets. In order to maintain the overall energy balance, 
the nuggets must emit energy at the same rate that it is 
deposited through proton annihilation, 



(1 - g)F ma = (l-g) 



dtdA ' 



(11) 



where 1 — g is the fraction of the annihilation energy that 
is thermalized. Note that both the rate of emission and 
the rate of annihilation are expressed per unit area A, 
so that the equilibrium condition is independent of the 
nugget size. The rate of annihilation F anu is 



F ann = 2 GeV • / ■ v ■ n VM (r) 



(12) 



Taking the typical values v ~ 10 c and 



300/ cm 3 



gives the relationship between the temperature and typ- 
ical parameters describing the nuggets, /, g, 



T 



4+1/4 



(1-.9) 



10" 



(14) 



As discussed in [4], reasonable values for / ~ 1/15 and 
g ~ 1/10 all lead toaT~ eV equilibrium temperature. 

The heat-capacity of the nuggets is estimated in Ap- 
pendix 5. If the gas of positrons occupies a substantial 
volume of the nugget, then the heat capacity is "large" in 
the sense that it will require many annihilations to raise 
the temperature of the nuggets to the eV scale at which 
equilibrium is established. Thus, the antimatter nuggets 
will act as effective thermal integrators, slowly reaching 
a relatively constant average temperature T ~ eV. 



where 2 GeV = 2 m p is the energy liberated by proton 
annihilation, v is the speed of the nugget through the 
visible matter, n VM (?) is the local visible matter density, 
and 



'arm -. „_ 1 

A 



is the factor by which the effective cross-section er ann for 
proton annihilation is reduced from the geometric cross- 
section A due to the possibility of reflection from the 
sharp quark-matter surface (in contrast, the positron dis- 
tribution in the electrosphere is very smooth), as dis- 
cussed in [4] . 

The typical galactic scale for the speed is v ~ 
100 km/s ~ 10~ 3 c, while the density at a distance r ~ 
kpc from the center is 



JT-VM 



f ■ n disk 



150 

cm 3 



where we have adopt a scaling behaviour close to that of 
an isothermal sphere [37] for the typical visible matter 



IV. EXPLAINING THE WMAP HAZE 

In our proposal, interstellar matter annihilates on an- 
timatter nuggets. The nuggets then radiate this energy 
over a wide range of frequencies. The model thus makes 
definite predictions relating these emissions: they should 
have similar morphologies, and the relative intensities 
should be related by an overall energy budget determined 
by the local annihilation rate. 

Four types of emission are from "hot spots" at the 
annihilation sites, and should be observable from the core 
of our Galaxy: 

(B.l) Electron annihilations through positronium pro- 
duce a well defined 511 keV emission [26, 27] 
that is consistent with, and could possibly explain 
the puzzling diffuse 511 keV emission observed by 

SPI / INTEGRAL. 

(B.2) Direct electron annihilation can also produce emis- 
sion in the 1-20 MeV band [28] which is consistent 
with, and could explain part of the diffuse gamma- 
ray emissions observed by COMPTEL. 
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(B.3) Proton annihilation produces keV x-ray emission 
from a hot spot at the annihilation site [4] that is 
consistent with, and could possibly explain the puz- 
zling diffuse x-ray emissions observed by Chandra. 
(B.4) Proton annihilation occasionally produces GeV 
photons [4] that are consistent with, and could 
partially account for the gamma-ray emissions ob- 
served by EGRET. 
All of these emissions are "direct" in the sense that the 
timescale for the emission is much shorter than the time 
between successive annihilations. Thus, the intensity of 
these emissions depends only 6 on the rate of annihilation 
events which is proportional to n VM (r)n OM (r) — the prod- 
uct of the local visible and dark-matter distributions at 
the annihilation site. The emitted spectrum is also in- 
dependent of the local density. We emphasize that the 
model makes two nontrivial predictions: 1) that the mor- 
phology of these emissions is very strongly correlated, 
and 2) that the spectral properties of these emissions are 
independent of position. 

A comparison between observations of the "direct" 
emissions B.l through B.4 along the same line-of-sight is 
possible because the local emission depends only on the 
local rate of annihilation <j>(r) oc n VM (r). The observed 
flux thus depends on the uncertain matter distribution 
through the same line-of-sight integral: 

$5H,x-ray,etc. OC J dQdl n VM (l)n DM (l) , (15) 

which cancels when comparing emissions from the same 
position in the sky. 

There is an additional emission from the 
nuggets: 

(B.5) Energy not directly released through one of the 
mechanisms B.1-B.4 heats the nuggets, ultimately 
being thermally radiated. 
As we have shown in Section III B and Appendix 5, 
the heat-capacity and energy budget ensure that the 
nuggets have a well-defined temperature scale of T ~ 
1 eV in the core of the Galaxy. The resulting thermal 
bremsstrahlung emission thus "averages" the annihila- 
tions over time, and the resulting emissivity and spec- 
trum will depend on the temperature T(n VM ) which is a 
function of the local visible matter density. 

The observed flux of this thermal emission exhibits a 
slightly different dependence because the local emission 
depends on the temperature </>(?) = <p{T[n VM (r)]} . As we 
shall show Appendix 6, however, the difference is small, 
and can be ignored for the order of magnitude compar- 
isons we present here. 

Finally, in principle, we may compare the total thermal 
emission (9) with the "direct" emissions because thermal 



equilibrium relates the rate of total emission to the rate of 
annihilation, both of which arc proportional to n VM (r). In 
practise, however, the thermal cV scale emission cannot 
be seen against the bright stellar background. 



A. Energy Budget 

In [4], the "direct" emissions B.l, B.3, and B.4 were 
compared, showing that our proposal is consistent with 
the current observations, and using the observations to 
constrain the properties of the nuggets. In particular, two 
parameters were introduced describe complicated proper- 
ties of the nuggets: The parameter / < 1 was introduced 
to describe the suppression of the proton annihilation 
rate with respect to the electron annihilation rate, and 
the parameter g < 1/2 was introduced to describe the 
fraction of the proton annihilation energy that is directly 
released as x-rays. We emphasize that these parameters 
are not free, but they depend on detailed models of the 
nuggets and are beyond the reach of present day calcu- 
lational techniques. 

By hypothesizing that emissions from the nuggets com- 
pletely explain the puzzling 511 keV (B.l) and diffuse x- 
ray (B.3) emissions, one obtains fg ~ 6 ■ 1CP 3 , which is 
consistent with the theoretical estimates, and provides a 
nontrivial test of the theory. 

Neglecting the small corrections to the line-of-sight av- 
eraging discussed above, we may perform a similar anal- 
ysis of the WMAP haze to see if it is also consistent with 
the our proposal. The thermal energy input into the an- 
timatter nuggets is the complementary fraction 7 1 — g 
of the total proton annihilation energy not directly re- 
leased as x-rays. Thus, if we use the observed x-ray flux 

Chandra to provide the energy normalization, then the 
total thermal emission will be approximately 



1-g 
9 



Chandra- 



(16) 



The total thermal emission $t may then be used to esti- 
mate the observed microwave emission in a specified fre- 
quency band by computing the ratio 7 of spectral emis- 
sivity (8) in the specified band to the total emissivity (9): 



7 = 



1 

Ftot 



hAo. dp 

dui — — (lo) 



25 - 121n(w/T cff ) 



60 T r 



Aw, 



cir 



where T e fT is an "average" temperature that accounts for 
variations along the line of sight. This is the fraction 
of the total emitted thermal radiation emitted in the mi- 
crowave band to of width Aw <C to <C T e g. For the typical 



There additional small dependencies that we neglect here: for 
example, on the local speed of the nuggets. This introduces 
only small uncertainties, however, and certainly do not affect 
the overall magnitude. 



7 Technically, we should include the energy from electron annihi- 
lations, and subtract the fraction a/a 3 of GeV photons B.4 [4]. 
These are only small corrections to the overall energy budget. 
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scale of the emissions we are considering, T e R ~ eV and 
co ~ h ■ 30 GHz ~ 10 -4 eV, so we have 



7 



T e ff 



The total observed microwave flux is then related to the 
total thermal flux (16) by 



Aw- 



d$WK 



1 



7" 
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Chandrai 



dui g 

giving us the relationship 

d$wMAP _ 2 
duj T eS g 

Observationally, Chandra observes a total flux [29] 

erg 



Chandra- 



Chandra ' 



2 x 10" 



(17) 



(18a) 



cur • s • sr 
while the WMAP haze flux is [5, 6, 7, 8] 

d$wMAP ,„ fi ,kJy (3-6) x lQ- 20 erg 

j = l-J-Dj 5 XJ • V l ° b ) 

auj sr cm z • s • sr • Hz 

Combining these, and converting 1 Hz « 4 x 10 -15 eV, 
we predict that the observed WMAP haze intensity will be 
saturated by thermal antinugget emission if the param- 
eters which enter in our estimate satisfy the following 
constraint 

eV l-g 
T c s g 

Although this relationship is only approximate, it is quite 
amazing that it is satisfied (without any adjustment) if 
the previous estimates for T e g and g are used. Thus, the 
nontrivial relationship (19) — which depends strongly on 
the observed intensity of the GHz WMAP haze — is satis- 
fied by the phenomenological parameters determined by 
considering only the keV scale emissions. 



(2-4). 



(19) 



B. Spectrum 

The observed spectrum of the wmap haze is extremely 
hard [5, 35]. This feature is easily accommodated in 
our model by the logarithmic dependence of the ther- 
mal bremsstrahlung emission (8). Indeed, the wmap 
haze was initially interpreted as thermal bremsstrahlung 
(free-free emission) from a hot (T ~ eV) gas [5] (10 4 K< 
T < 10 6 K), but this interpretation was rejected because 
a Ha recombination line, which should accompany the 
haze, is not seen. (The possibility of much hotter plasmas 
T > 10 4 K has also been ruled out [5, 6, 7, 35].) 

It is quite remarkable that the T ~ 10 4 eV scale 
arises natural in our proposal in two completely indepen- 
dent ways: (14) and (19). Our proposal, thus, naturally 
fits the observed spectrum of the WMAP haze, without 
any Ha recombination line since the emission is from a 
purely positronic gas. It is also remarkable that the spec- 
trum exactly corresponds with bremsstrahlung radiation, 
as was originally suggested in [5] to fit the data. 



C. Morphology 

Our proposal also makes a definite prediction about the 
morphologies of the various emissions. In particular, the 
morphology of the "direct" emissions B. 1-B. 4 should 
be almost identical. As discussed in Appendix 6, even 
though the wmap haze is a thermal emission, the de- 
pendence on the line-of-sight matter distribution is quite 
weak, and our model thus also predicts that the mor- 
phology of the WMAP haze be closely correlated with the 
morphology of the "direct" emissions. 

Ultimately, if our proposal is correct, the morphology 
of all of the emissions are direct probes of the matter dis- 
tributions in our Galaxy, and may thus become a useful 
measuring tool. 



V. SUMMARY 

We have now demonstrated that our proposal naturally 
and nontrivially explains diffuse galactic radiations B.l- 
B. 5 over thirteen orders of magnitude from microwave 
(10~ 4 eV) to GeV scales. The only "parameter" in our 
model is the overall size of the nuggets, and the depen- 
dence of the observations on this parameter can virtu- 
ally be eliminated by comparing observations along the 
similar lines-of-sight to the core of the Galaxy. Com- 
parisons along the same line-of-sight also virtually elimi- 
nate uncertainties related to the distribution of matter in 
our Galaxy. Removing this uncertainty, our proposal de- 
pends on only a couple of parameters: / < 1 and g < 1/2. 
These represent presently intractable calculations, but 
have tight upper bounds and can vary by no more than 
an order of magnitude or so. 8 

Together with our previous work [4], we now have 
several constraints on these parameters by postulating 
that matter annihilation on dark-antimatter nuggets ex- 
plains significant proportions of puzzling diffuse emis- 
sions. These constraints provide a highly nontrivial test 
of our proposal. 

(S.l) The diffuse 511 keV emission B.l observed by 
SPi/ integral has been identified as primarily due 
to positronium annihilation. The puzzle is how 
positrons come to be diffusely distributed through 
the core of our Galaxy. We propose that the 



8 The parameters / and g can, in principle, be calculated from the 
first principles. However, such a computation is extremely diffi- 
cult as it requires solving the many-body problem in a strongly 
coupled regime. Indeed, the phase diagram of quark matter in 
the relevant regime is still largely unknown. Even with these 
reservations, we still are quite confident that these parameters 
cannot deviate much from their "natural" values. In this respect 
our proposal is predictive: there is little freedom to change these 
parameters. This is in contrast with most other dark matter pro- 
posals where parameters can be arbitrary changed by many or- 
ders of magnitude to exploring an enormously large, and largely 
unconstrained parameter space. 
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positron clcctrosphere of dark-antimatter nuggets 
provide the source of positrons. This source is dis- 
tributed diffusely, and produces a definite spectrum 
consistent with the observed spectrum [26, 27] that 
is independent of any model-specific parameters. 
One prediction is that the spectrum is independent 
of observed direction. Another is that the intensity 
is determined by the product n VM (r)n DM (r) of the 
distribution visible and dark matter. We use this 
emission as a baseline to which we compare the 
other emissions in order to remove the uncertain- 
ties of the nugget size and the matter distribution 
along the line of sight to the galactic core. 

(5.2) Associated with this emission is the MeV spectrum 
from direct c + c~ annihilation B.2. This spectrum 
is model independent, and consistent with observa- 
tions and background models, possibly explaining 
the 1-20 MeV energy deficits seen in the COMPTEL 
data [28]. 

(5.3) The diffuse keV x-ray emissions B.3 measured by 
Chandra arc puzzling because of the implied energy 
budget. The spectrum looks like a thermal 8 keV 
plasma [29], but such a plasma would not even 
be gravitationally bound and would require a huge 
unidentified source of energy to fuel. We propose 
that this emission is due to bremsstrahlung emis- 
sion from positrons excited from protons annihilat- 
ing on the dark-antimatter nuggets [4]. The spec- 
trum for this process is also largely independent 
of model-specific parameters, is consistent with the 
observations, and is also independent of position. 
Comparing this emission the 511 keV emission gives 
one constraint on the parameters fg [4], 



fg w 6 x 10- 



(20) 



that is satisfied by their natural scales / f=a 1/15 
and g ~ 1/10. The morphologies should also be 
related, depending on the product n VM (r)n DM (r). 

(5.4) The direct emission of GeV photons from proton 
annihilation B.4 is consistent with gamma ray ob- 
servations by egret, explaining up to one tenth or 
so of the observed spectrum [4] . 

(5.5) The annihilation energy not released through one 
of the previous sources of emission ultimately ther- 
malizes. As shown in Section IIIB, the nuggets 
reach equilibrium with a typical T ~ cV scale 
through the constraint (14), 



T 
vV 



4+1/4 



/ 



io- 



(14) 



which is again satisfied by the natural parame- 
ter scales. We emphasize that this constraint is 
virtually model independent, depending on only 
the emissivity calculated in Section III A and the 
properties of the matter distribution in the core of 
the Galaxy. The emissivity is dominated by the 
properties of the nugget electrosphere in the low- 



density regime, where calculations are under ordcr- 
of-magnitude control. Although not particularly 
sensitive to the emissivity, the scale set by (14) 
would be altered by an order of magnitude if the 
emissivity were black-body. 



A. Predictions 

The constraints described in Section V did not involve 
any measurements of the wmap haze. Thus, the temper- 
ature T ~ eV (14) allows us to unambiguously predict 
the energy emitted in the tail of the thermal distribution 
in the microwave band. Amazingly, comparing this with 
the observed WMAP haze gives (19) 



eV 1 

Teg 
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(2-4), 



(19) 



which is satisfied with the natural values of the param- 
eters, even though the frequencies are at a scale many 
orders of magnitude smaller than the scale at which the 
parameters were constrained. 

Unlike (14), this estimate is extremely sensitive to the 
flat spectral properties of the thermal emission (8) at 
small frequencies u) <C T: a very specific feature of ther- 
mal bremsstrahlung emission. If the spectrum were not 
flat, there is no way that this constraint would be satis- 
fied. 

Another difference between the two estimates (14) and 
(19) is that the former is sensitive overall normalization of 
(9) whereas the latter is sensitive to the detailed shape 
of the spectrum (8). We also emphasize that the con- 
straints (14) and (19) not only deal with emissions from 
scales separated by 8 orders of magnitude, but that con- 
dition (19) depends on observed intensity of the wmap 
haze which has absolutely nothing to do with the esti- 
mate (14). That these two estimates agree with each 
other and with the estimates in [4], (where the intensity 
of the diffuse keV x-rays were compared with the 511 keV 
emission), is truly remarkable. 

In addition to satisfying the constraints described in 
the previous section, our proposal also makes the defi- 
nite prediction that the morphologies of the 511 keV flux, 
the 1-20 MeV 7-ray emission, and the x-ray flux should 
all be identical, following the distribution n VM (f)n DU (r) 
of visible and dark matter. The morphology of the 
wmap should be similar but may differ slightly because 
of the different line-of-sight integral (40). For example, 
Chandra has detected a diffuse x-ray emission with flux 
6.5 x 10" 11 erg/cm 2 /s/dcg 2 from a region of the disk 28° 
off the center [38] . This is one order of magnitude smaller 
than the observations from the core of the Galaxy, and 
so our model predicts that the microwave emission from 
this region is about one order of magnitude smaller than 
the WMAP haze from the galactic core, which seems con- 
sistent with the observations [5]. 

Finally, this proposal makes definite testable predic- 
tions for the properties of the emitted spectra. In par- 
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ticular, all bands but the WMAP haze are produced on 
an event-by-event basis, and are thus independent of the 
rate at which the annihilation processes occur. The ob- 
served spectra should thus be largely independent of the 
direction of observation. Only the intensity should vary 
as a function of the collision rate, and this should be cor- 
related with the visiblc/dark-matter distribution as dis- 
cussed above. The WMAP haze will have a slight spectral 
dependence, but only through the temperature depen- 
dence T(n VM ) oc fivd 17 which is quite weak. 

B. Conclusion 

Our dark matter proposal not only explains many as- 
trophysical and cosmological puzzles, but makes definite 
predictions about the correlations of the dark and visi- 
ble matter distributions n VM Tiim with five different bands 
of radiation ranging over 13 orders of magnitude in fre- 
quency. In addition, it makes the definite prediction that 



these spectra of the emissions should be virtually inde- 
pendent of the local environment. Such correlations and 
spectral properties would be very difficult to account for 
with other dark-matter candidates. Future observations 
may thus easily confirm or rule out this theory. If con- 
firmed, it would provide a key for many cosmological and 
astrophysical secrets, and finally unlock the nature of 
dark matter. 
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APPENDIX: DETAILED CALCULATIONS 
1. Structure of Electrosphere 

Here wc briefly discuss the properties of the antimatter 
nugget electrosphere to set the scales of the problem and 
determine the structure required in Section III A. 

The radius of the nuggets depends on the mass, but 
must be larger than R > 10 cm at the lower limit 
|£>| > 10 20 set by terrestrial nondetection. We ex- 
pect that the most likely size is on the order of \B\ ~ 
10 25 " 27 [18]. The quark-matter core of the nuggets ends 
sharply on a fm scale as set by nuclear physics. Near 
the surface, as the density falls, the quark matter will 
definitely be charged due to the relatively large mass 
of the strange quark m s ~ 100 MeV [39, 40, 41, 42], 
however, depending on the phase of quark matter real- 



ized in the core, the matter may be charged through- 
out. Charge neutrality will be maintained through beta- 
equilibrium, which will establish a positron chemical po- 
tential fi c + = /xrj — 10 MeV. (The precise value de- 
pends on specific details of the quark-matter phase and 
may range from a few MeV to hundreds of MeV, but is 
about an order of magnitude less than the quark chem- 
ical potential [i q ~ 500 MeV [39, 43].) This will induce 
a thin but macroscopic "electrosphere" of positrons sur- 
rounding the quark-matter core in the transition region 
as fx e + — > in the vacuum. 

The structure of this electrosphere has been consid- 
ered for quark matter [39, 40], and the existence of this 
"transition region" is a very generic feature of these sys- 
tems. It is the direct consequence of Maxwell's equations 
and chemical equilibrium. The region is called the elec- 
trosphere, emphasizing the fact that quarks and other 
strongly interacting particles arc not present. In the 
case of antimatter nuggets the "electrosphere" comprises 
positrons. 

The variation of this chemical potential (i Q +(z), and 
the density n(z) as a function of distance z from the 
surface of the nugget may be computed using a mean- 
field treatment of the Maxwell equations [39, 40, 44]. For 
example, in the relativistic regime, one has [42] 



H c +{z) 



2a (z + z ) ' 



n[Z) ~ 3^ 2 ~ 3^ 2 I 2a 



(z + zo) 3 ' 



(21) 




where (j,q = n c +(z = 0) ~ 10 MeV is the chemical po- 
tential realized in the nugget's bulk. The corresponding 
results can be obtained outside of the relativistic regime, 
but they do not have a simple closed form. These calcu- 
lations treat the electrosphere as a one-dimensional wall 
rather than including the full radial structure, essentially 
keeping only the first term in the z/R expansion. This 
approximation does not affect the order of magnitudes of 
our calculation and will also be employed here. 

The majority of the thermal emission considered in 
this work comes from the nonrelativistic regime, which 
we may also analyze analytically using the Boltzmann 
approximation. The mean-field approximation amounts 
to solving the Poisson equation 



-47ren(?) 



(22) 



where 0(r) is the electrostatic potential and n(r) is the 
density of positrons. Using the spherical symmetry of the 
nuggets and making the one-dimensional approximation, 
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we can write this as 9 

d 2 0(z) 
dz 2 



-47ren(z) 



(23) 



where z is the distance from the quark nugget sur- 
face. We now introduce the positron chemical poten- 
tial n e +{z) = —e(j)(z) which is the potential energy of a 
charge at position z relative to z = oo where we take 
H e + (oo) = as a boundary condition. This gives 



d 2 fx e+ (z) 



Awa n[^ c + (2)] 



(24) 



regime: 



n{z = 0) 



T 



2naz 2 



(mT) 



3/2 



(28) 



While not an exact matching procedure, this will give the 
correct parametric dependence and will be valid for the 
order-of-magnitude estimates required. Thus, we have 
the following characterization of the density in the Boltz- 
mann regime as used in Section III A: 



with the additional boundary conditions /j, c + (z = 0) = 
Ho ~ 10 MeV as established by beta-equilibrium in the 
quark matter. Here n[/i G +] is the density of a free Fermi- 
gas of positrons as a function of the chemical potential 
for the positrons. The full relativistic form is 



n[n] 



where r. 



d 3 p 



1 



1 



1 1 ""^ 
1 + e t 



1 + e t 



(25) 



which has the property n[/i 



0] = required by our identifying the chemical potential 
/i with respect to z = 00. This is quite complicated, but 
is well approximated in the nonrelativistic Boltzmann 
regime where n <C (mT) 3 / 2 by 



(2^ 



^(n^y'^/T (26) 



where we have performed a nonrelativistic expansion 



tn+p 1 /(2m), dropped the antiparticle con- 



tribution, and neglected the quantum degeneracy. The 
effective chemical potential fi = /i c + — m is related to the 
vacuum chemical potential \i by subtracting the mass. 
We note that the right boundary condition must now 
be changed to n(z = 00) = because jl does not tend 
to zero under these approximations. The left boundary 
condition must be determined by matching the density 
at some point to the full relativistic solution that inte- 
grates to the quark-matter core and matches (21). The 
differential equation (24) may now be expressed in terms 
of either p>(z) or n(z), leading to the peculiar solution 



n(z) 



T 



1 



2na (z + z) 2 



(27) 



where z is an integration constant determined by match- 
ing to a full solution. Here we make a simple approx- 
imation, defining z = as the onset of the Boltzmann 



9 Here we drop the radial term 2<f>'(z)/r on the left hand side of 
(22), assuming that the radius of the nuggets K > 2 is much 
larger than the thickness of the electrosphcre. 



z 1 ~ \[2~tux ■ m ■ \ — . 

V m 



(29) 



The region where z < here corresponds to the region 
of higher density closer to the nugget's surface where the 
Boltzmann approximation breaks down due to degener- 
acy effects. One can argue, however, that in this degener- 
ate regime, the emissivity is strongly suppressed for two 
reasons: 1) only a small portion of the particles close to 
the Fermi surface can participate in scattering processes, 
so the phase space for emission is dimensionally reduced, 
and 2) as the density increases, the plasma frequency 
increases, and low-energy photons cannot escape. 

Our estimates include only emission from the Boltz- 
mann regime. In principle, the emission from denser re- 
gions could contribute at the same order: this approxi- 
mation thus underestimates the emission by a factor, but 
not by an order of magnitude. 

We make one final set of remark in response to 
the criticism [45] (version 3) of our proposal. As we 
have shown, the density profile behaves very differently 
in different physical regimes. In the ultrarelativistic 
regime (21), one has the dependence n oc V 3 where V(z) 
is the electrostatic potential in the mean-field approxima- 
tion [39, 40, 44], whereas in the nonrelativistic Boltzmann 
regime (26), one has n oc e v ^ T , which is the well-known 
expression for the density in nonrelativistic and nondc- 
generate systems [46]. In the intermediate regimes, the 
dependence is quite complicated due to the competing 
scales that appear in n[n] (25). Thus, one cannot sim- 
ply apply formulae like n (x V 3 derived in one regime to 
describe physics in another as was done in [45]. 

In general, one must also take into account quantum 
many-body effects — such as charge screening (completely 
ignored in [45]), the plasma frequency, etc. In the Boltz- 
mann regime discussed here, the density is sufficiently low 
that many-body corrections may be neglected and vac- 
uum results, such as the cross-section (2), employed. At 
higher densities, however, when the degeneracy becomes 
important — of order roug hly n > (mT) 3 / 2 m our case — 
many-body effects can drastically alter the behaviour of 
the system and cannot be neglected as they were in [45] , 
even when considering only the qualitative physics. 
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2. Boltzmann Averages 

To evaluate equation (3) we need to perform the ther- 
mal average 



vi -V 2 \ 17+ 12 In 



(Pi - P2) 



(30) 



As we are in the Boltzmann regime, we may simplify the 
calculation by computing this in the Boltzmann ensem- 
ble. Formally, we must integrate over both momenta p\ 
and P2 , but as we are only interested in the order of mag- 
nitude, we simply perform the average over only a single 
momentum p 2 , setting pi = 0. 
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where 



h{x) = 17+ 12g(x), 
g(x) = ln(2) + 

Ei(x) -- 



1 + x 
-dz. 



The following approximation for h(x) is accurate to 
within 25% for all x: 



h(x) 



17-121n(x/2) x < 1, 
17+121n(2) x>l. 



(32) 



3. Finite-Size Effects 

So far, our calculations have assumed that we are work- 
ing in infinite matter. Here we estimate the size of the 
corrections due to the fact that the nuggets have a finite 
extent on the order of L > 10~ 5 cm. We shall demon- 
strate that properly accounting for these corrections does 
not significantly affect our estimates of the microwave 
band emission (though it may drastically suppress emis- 
sion at much longer wave lengths). 

In principle, finite-size effects may change the cross sec- 
tion (2), and therefore, our estimation of the emissivity 



10 Including the full angular integrals increases the result by a factor 
of about v2 or so. 



vi2 17+ 12 In ■ 



(^) 



'/ T l 



17- 



mir \ TV V " \TJ 

The calculation of g(x), however, is somewhat tricky. 



(3). The cross-section (2) was derived using a continuum 
of plane-wave states, whereas to account for the finite- 
size effects, one should use the basis of states bound to 
the quark core. To estimate the size of the corrections, 
one can imagine confining the positrons to a box of finite 
extent L. 

The electromagnetic field may still be quantized as in 
free-space with states of arbitrarily large size because the 
photons are not bound to the core, and are not in ther- 
modynamic equilibrium with the positrons. Their mean- 
free-path is much larger than L, so the low-energy pho- 
tons produced by the mechanism described above will 
simply leave the system before they have a chance to 
interact with other positrons. 

Therefore, it is only the positron states that must be 
considered over a finite-size basis, which will modify the 
corresponding Green's function used in the calculation of 
the cross-section (2). These modifications occur for mo- 
menta of the scale 5p ~ H/L. If L > 1CP 5 cm, then this 
corresponds to shifts in the energies of SE ~ (Sp) 2 /2m ~ 
10~ 6 cV -C 10~ 4 eV, which is much smaller than the 
transitions responsible for the emission at microwave fre- 
quencies. Thus, we conclude that finite-size effects do 
not drastically change the positron Green's function in 
the region of interests. In other-words, the expression 
for the cross-section (2) — derived using the standard (in- 
finite volume) Green's functions — remains valid for the 
estimation of the emission and spectrum down to the 
microwave region 10 -4 eV. We also note that finite-size 
effects do not change our estimates for the density (25, 
26) because the finite-size effects SE <C T are much 
smaller than the typical energetic scale T ~ cV of the 
problem. Thus, our expression (2) remains valid for pho- 
ton energies u) > 10~ 4 eV. To calculate the emission 
of radiation with much longer wavelengths, however, re- 
quires one to account for these finite-size corrections, and 
we expect the emission of extremely low-energy photons 
<jj <C SE < 10~ 6 cV to be suppressed. 

One may ask how microwave radiation may be emitted 
from the nuggets when the wavelength A is much larger 
than the size of the nugget A L. 11 In general this 
is not a problem — consider the well-known astrophysical 
emission of the A = 21 cm line from hydrogen with a size 
a ~ 10 -8 cm — but there is a potential suppression. The 
coherence time t of the positrons which must be com- 
pared with the formation time ~ lu~ 1 of the photons. If 
the coherence time is too short, then multiple scatter- 
ings will disrupt the formation of the photons. This sup- 
pression is a case of the so-called Landau-Pomeranchuk- 
Migdal (lpm) effect. 

To estimate the coherence time t, consider the cross- 
section (Tee of the positron-positron interaction. This 
scales as <j ee ~ a 2 /q 2 where q ~ 6 _1 is the typical 
momentum transfer, and may be expressed in terms of 



We thank a referee for raising this question. 
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the impact parameter b ~ n -1 / 3 , which is estimated in 
terms of average interparticle spacing where n is the local 
positron density. 



The mean-free-path I is thus I 1 ~ cr ee n 



2„l/3 



an 



Therefore, the typical time between collisions (which is 
the same as coherence time) is r ~ l/v where v ~ yjT/m 
is the typical positron velocity. 

Collecting all factors together we arrive at the estimate 



2 nV3 V T a 2 T 



> 1. 



(33) 



It is clear that this condition is satisfied for uj > 10 -4 eV 
and T < 1 eV, even for z = 0. Thus, we were justi- 
fied in omitting LPM effect in our estimates in the low- 
density regime (27). However, from the same estimate 
it is clear that this suppression becomes important for 
either smaller frequencies uj <C 10~ 4 eV or at higher den- 
sities. We shall now show that there us a much more 
significant high-density suppression than the lpm effect, 
which effectively turns off emission from the bulk of the 
nuggets where the positron density is significantly higher 
than (27). 



including these effects: 

ni(z,T)n 2 (z,T) -> 
d 3 pi d 3 p 2 6(i 



Uj)8(uj — LU p ) 



(2tt) 3 (2tt) 2 



1 + e 



e p2 



(34) 



where e Pl and e P2 are the colliding positron quasipar- 
ticle energies above the Fermi surface. The factor of 
9(e Pl + £ P2 — uj) accounts for energy conservation: the 
initial energy must be larger than the energy of emitted 
photon; and the factor 9(u> — u) p ) accounts for the effects 
of the plasma frequency. Only photons with uj > uj p can 
propagate in the dense media: photons with smaller en- 
ergies will be absorbed on distances ~ uj" 1 . As we shall 
see, this leads to an exponential suppression of emission 
~ exp(— ujp/T) when uj p 3> T. As such, we have omitted 
the aforementioned lpm suppression and additional Pauli 
blocking effects, as these are comparatively insignificant. 
The integral in (34) can be estimated from 



d£p 1 d6p 2 



T 2 e 



4. Emission from Very Dense Regions 

Here we estimate the emissivity from very dense re- 
gions of the nugget when the Boltzmann regime breaks 
down and degeneracy plays a crucial role. As we shall 
argue below, the corresponding emission from very dense 
regions can be neglected in comparison with the esti- 
mates (8) and (9) used in the text. 

We start from (3), which remains valid for any densi- 
ties n(z,T). To deal with denser regions, however, one 
must the full expression for (25) which includes the ef- 
fects of quantum degeneracy. In these dense regions, 
only the states close to the Fermi surface are excited and 
can participate in emission: the states deep within the 
Fermi surface are "Pauli blocked" and cannot partici- 
pate in low-energy interactions. It is the density of these 
"quasiparticles" — not the full density — which enters the 
emissivity calculations. The other key property for these 
estimates is the is the plasma frequency uj p , which char- 
acterizes the propagation of photons in the degenerate 

systems. For ultrarelativistic systems, uj p = -ff- while 
for nonrelativistic systems, uj p — 4 ™ n . The plasma fre- 
quency can be thought as an effective mass for the pho- 
ton: only photons with energy larger than this mass can 
propagate outside of the system. Photons with u> < uj p 
are "off-shell" or "virtual" : these can only propagate for 
a short period of time/distance ~ uj~ 1 before they decay. 

To derive the analogues of equations (8) and (9) for 
the denser regions, we must start with (3), but insert 
the proper form for the expression for n\(z,T)n2{z,T), 



to give 

n l{z ,T)n 2 (z,T)^ P ^^e^(-^). (35) 

The main point is that tu p grows with the density as 
Lj p ~ y/n in the nonrelativistic regimes, and as uj p ~ <fn 
in the relativistic regimes. This leads to an exponen- 
tial suppression of emission from the dense regions of 
the nugget. The suppression is lifted only when uj p ~ T 
which occurs (as can be verified numerically) only when 
the densities are sufficiently low that the Boltzmann ap- 
proximation is valid. Our estimates (8) and (9), which 
are based on this approximation, are thus justified. 

One can estimate that the plasma frequency uj p is a 
few eV for densities (27) typical of the Boltzmann regime. 
Given our previous discussion, one might ask: How can 
low-energy photons lu < lo p still be emitted? The rea- 
son is that, although these photons would be reabsorbed 
in infinite matter, this reabsorption happens on a length 
scale of ujp 1 . At the typical densities in the Boltzmann 
regime, lo^ 1 ~ 0.3 ■ 10 -5 cm is sufficiently large com- 
pared with the nugget size that many of these photons 
will have left the nugget before being reabsorbed. One 
can interpret this effect as a decay of a quasistationary 
state in quantum mechanics through the tunnelling pro- 
cess, where the barrier has a size comparable with inverse 
energy of the quasistationary state. The suppression fac- 
tor only becomes sufficiently large when cu p > L" 1 . Then 
the reabsorption happens well before the photons emerge. 
This is the origin of the primary suppression in the very 
dense regions we have just discussed. 
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5. Heat Capacity 

Here we make a rough estimate of the heat capacity 
of the nuggets. At eV temperatures, the only modes 
that can be excited are the neutral Nambu-Goldstonc 
supcrfluid mode, which contributes as T 3 , and the gapless 
positrons at their Fermi surface of pp f=a 10 MeV. It is the 
gapless positrons that will dominate the low-temperature 
heat capacity of the nuggets. 

To estimate the heat capacity, we simply count the 
number of low-energy modes. The heat capacity density 
for a single mode with dispersion e p is 



c v 



1 



d 3 p d 

(2?r)3 ep JT \e^' T ± 1 



(36) 



The total heat capacity Cy ~ Vcy is obtained by inte- 
grating this over the volume V . Here are some approxi- 
mate contributions to the heat capacity from the nuclear 
core, assuming the density p ~ 5 GeV/fm 3 is a few times 
nuclear density so that V f» 26 B GeV -3 , and taking 
T w 1 cV: 

Single Boson: For a single boson, e p = p/3. This gives 
cy w 35 T 3 and a total contribution of Cy ~ 
10" 24 B. 

Two Fermions: For a pair of positrons with Fermi sur- 
face pf we have e p ~ v \p — pp | which gives cy ~ 
PpT/(3v) and a total contribution of Cy ~ 10 -12 _B 
where we take the Fermi velocity v w c for relativis- 
tic systems and ~ 10 MeV. 

In the complete absence of positrons, the heat capacity of 
a medium nugget B ss 10 24 would be about unity, mean- 
ing that the addition of 1 GeV of thermal energy would 
raise the temperature by a GeV, however, even if the 
core of the nuggets is completely neutral (for example, 
if the Colour-Flavour-Locked (cfl) phase were realized 
with equal numbers of up, down, and strange quarks), 
the surface of the quark-matter core will still be charged, 
requiring a large number of positrons in the electrosphere 
to neutralize the objects. 

Once one includes the contribution of the positrons, 
however, the heat capacity becomes much larger, and 
the input of 1 GeV of energy will only effect a very 
small change in temperature. For example, if positrons 
are present in the core, even the very smallest possible 
nuggets 12 B ~ 10 21 would heat by only an eV given 
1 GeV of energy We expect the nuggets have a typi- 
cal size of B io 25 " 27 [18], so nuggets would require 
many collisions in order to reach the eV temperature 
scales at which they would thermally equilibrate by bal- 
ancing the annihilation energy input with the thermal 
radiation (14). 



12 The constraint that the average B > 10 21 is quite strong and 
comes from negative terrestrial based search results [47]. 



Note, however, that while the heat capacity is large in 
the sense that many collisions are required to heat the 
nuggets to the equilibrium cV scale, it is small enough to 
allow the nuggets to heat to this temperature quickly on 
galactic scales. To see this, consider the total annihilation 
rate (12). 



1 aim i , - 1 1 aim 

at 



\ 10 25 J 



2/3 



GeV 



(37) 



where A is the surface area of the nuggets 

/3BGcV\ 2/3 ,/,„ o 

A = 4tt(-^——\ = 1.65 J B 2/3 fm 2 . (38) 

Therefore, the typical time between annihilations is on 
the order of 1 per second for nuggets with B ~ 10 25 . 
A nugget with positrons throughout (corresponding to a 
larger heat capacity) would require 



the 



1 eV • C v 



B \ 1/3 
— I hours (39) 



to reach eV temperatures. Thus the nuggets reach their 
equilibrium quite rapidly once they enter a region of high 
visible matter density. 



6. Line-of-Sight Integration 

A comparison between the "direct" emissions B.1-B.4 
is facilitated by the fact that they depend on the same 
line-of-sight average (15) 

$5H,x-ray,etc. °C / dfidZ Tl VM (l)n DM (l). 

In principle, comparing these with thermal emissions is 
more difficult because the line-of-sight integral has an ad- 
ditional dependence on the visible matter distribution. 
This arises because the emissivity depends on the tem- 
perature T of the nuggets, which in turn depends on the 
annihilation rate through the thermal equilibrium condi- 
tion (14). 

As can be seen from (9), (11) and (12), the temperature 
of the nuggets is related to the local visible matter density 
as T(n VM ) oc n™ 17 , whereas from (8), we see that the 
microwave emission scales as <fi(T) oc 7 113 / 4 oc n*, 3 / 17 = 

?1vm"-vm^ 17 - Thus, the microwave emission depends on 
the line-of-sight integral 



$wmap oc / dQdl [n VM (0«D M (0] n vM /17 (0- 



(40) 



Of course, one must also account for rescattering, ab- 
sorption and other effects along the line-of-sight, but our 
point is that the difference between the line-of-sight aver- 
aging of the "direct" emissions and the "thermal" emis- 
sions is highly suppressed — depending only on n v ^/ 17 . 
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This allows us to directly compare all emissions in order 
to estimate the order of magnitude for the energy budget 
discussed in Section IV A. Once the other calculations 
and observations are better constrained, one might be 
able to search for this type of scaling to test our proposal 
more rigorously. 



